Device and method for measuring phase retardation distribution and fast axis azimuth angle distribution in real time

ABSTRACT

Device and method for measuring phase retardation distribution and fast axis azimuth angle distribution of birefringence sample in real time. The device consists of a collimating light source, a circular polarizer, a diffractive beam-splitting component, a quarter-wave plate, an analyzer array, a charge coupled device (CCD) image sensor and a computer with an image acquisition card. The method can measure the phase retardation distribution and the fast axis azimuth angle distribution of the birefringence sample in real time and has large measurement range. The measurement result is immune to the light-intensity fluctuation of the light source.

CROSS-REFERENCE AND RELATED APPLICATIONS

The subject application is a continuation-in-part of PCT international application PCT/CN2012/000990 filed on Jul. 23, 2012, which in turn claims priority on Chinese patent application No. CN 201210199435.4 filed on Jun. 15, 2012. The contents and subject matter of the PCT and Chinese priority applications are incorporated herein by reference.

TECHNICAL FIELD OF THE INVENTION

The subject application relates to polarization measurement, especially one device and method for measuring phase retardation distribution and fast axis azimuth angle distribution of the birefringent component in real time.

BACKGROUND OF THE INVENTION

Birefringent component is widely used in many domains, such as polarized illumination system of immersion lithography, phase shifting interferometry and biological optics. The phase retardation and fast axis azimuth angle are two important parameters of birefringent component. The phase retardation distribution and fast axis azimuth angle distribution of the birefringent component must be acquired when it is used in the immersion lithography polarized illumination system and phase shifting interferometry. So it is of essentiality to precisely measure the phase retardation distribution and fast axis azimuth angle distribution of the birefringent component.

The Chinese Laid-Open Patent Application No. 200710178950.3 discloses a method and system of precisely measuring the optical phase retardation. By using an optical modulator in the optical setup, the modulated polarized light is generated. After filtering the measurement signals, the measurement of DC zero value is transferred to the measurement of AC zero value. By precisely judging the position of the extreme point, the phase retardation can be measured. But this method cannot measure the fast axis azimuth angle of the sample and the phase retardation distribution in real time.

Tsung-Chih Yu et al. in a paper entitled “Full-field and full-range sequential measurement of the slow axis angle and phase retardation of linear birefringent materials”, Applied Optics, Vol. 48, p. 4568(2009), discloses a method of measuring the phase retardation distribution and fast axis azimuth angle distribution of birefringent materials by using heterodyne interferometry method and three-step time-domain phase shifting method. This method needs to change part of the light path to measure the phase retardation distribution and fast axis azimuth angle distribution in step by step and uses the time-domain phase shifting technology, thus it is not feasible for measuring the phase retardation distribution and fast axis azimuth angle distribution in real time.

SUMMARY OF THE INVENTION

The purpose of the subject invention is to overcome the shortages of the above technology. One method and device of measuring the phase retardation distribution and fast axis azimuth angle distribution in real time is proposed. The measured result is immune to the fluctuation of the initial light intensity and this method has large measuring range.

The subject application provides a device for measuring the phase retardation distribution and fast axis azimuth angle distribution in real time, which is composed of a collimating light source, a circular polarizer, a diffractive component, a quarter-wave plate, an analyzer array, a CCD image sensor and a computer equipped with an image acquisition card. Said analyzer array is composed of four analyzers whose polarization directions successively increase by 45°. They are respectively named as the first analyzer, the second analyzer, the third analyzer and the fourth analyzer. The positions of the above components are as follows:

Said quarter-wave plate is located in the same path with said first analyzer. And the angle between the fast axis of said quarter-wave plate and the transmission direction of said first analyzer is 45° or 135°. Light emitted from said collimating light source passes through said circular polarizer and said diffractive component and then is split into four sub-beams. One sub-beam passes through said quarter-wave plate and is then analyzed by said first analyzer. The other three sub-beams are directly analyzed by said second analyzer, said third analyzer and said fourth analyzer, respectively. The output port of said image sensor is connected to the input port of said computer. The faucet for the measuring sample is set between said circular polarizer and said diffractive component.

Said collimating light source is a He—Ne laser, solid state laser or semiconductor laser.

Said circular polarizer is made of birefringent crystal, birefringent film or micro optical element.

Said diffractive component is a Quadrature Amplitude grating, Quadrature Phase grating or Dammann grating. The diffractive component can split the incident beam into four sub beams with the same light intensity.

Said quarter wave plate is a crystal wave plate, prismatic wave plate, film wave plate or composite wave plate.

Said first analyzer, second analyzer, third analyzer and fourth analyzer are all the polarizer whose extinction ratio is better than 10⁻³.

The method of measuring the phase retardation distribution and fast axis azimuth angle distribution in real time contains the following steps.

{circle around (1)} Insert the measuring sample into said faucet between said circular polarizer and said diffractive component, then adjust the light beam to let it perpendicularly pass through the measuring sample.

{circle around (2)} Turn on said collimating light source, said CCD image sensor and said computer. Said CCD image sensor receives an image formed by said four sub-beams and then transmits it to said computer. Said computer segregates said image into four sub-images. Then the four said sub-images are all pixelated and established a same coordinate system, respectively. The measuring sample is also matrixed and established a coordinate system which is the same with that of the sub-image. The intensity values in said four sub-images corresponding to the measured matrix unit (x, y) of the measuring sample are I₁(x, y), I₂(x, y), I₃(x, y) and I₄(x, y). By calculating said intensity values, the phase retardation and fast axis azimuth angle of said matrix unit (x, y) of the measuring sample are obtained. So the phase retardation distribution and the fast axis azimuth angle distribution of the measuring sample can be obtained by processing said four sub-images.

Said computer processing said sub-images contains the following steps in detail: When the polarization direction of said first analyzer is 45° relative to the fast axis of said quarter wave plate, said computer will perform the steps {circle around (3)}, {circle around (4)}.

{circle around (3)} Said computer processes the intensity values I₁(x, y), I₂(x, y), I₃(x, y) and I₄(x, y) corresponding to said matrix unit (x, y) as follows:

${V_{1}\left( {x,y} \right)} = {{{\sin \left( {\delta \left( {x,y} \right)} \right)}{\sin \left( {2{\theta \left( {x,y} \right)}} \right)}} = {\frac{2\; {I_{1}\left( {x,y} \right)}}{{I_{1}\left( {x,y} \right)} + {I_{3}\left( {x,y} \right)}} - 1}}$ ${{V_{2}\left( {x,y} \right)} = {{{\sin \left( {\delta \left( {x,y} \right)} \right)}{\cos \left( {2{\theta \left( {x,y} \right)}} \right)}} = {\frac{2\; {I_{4}\left( {x,y} \right)}}{{I_{1}\left( {x,y} \right)} + {I_{3}\left( {x,y} \right)}} - 1}}},{{V_{3}\left( {x,y} \right)} = {{\cos \left( {\delta \left( {x,y} \right)} \right)} = {1 - {\frac{2\; {I_{2}\left( {x,y} \right)}}{{I_{1}\left( {x,y} \right)} + {I_{3}\left( {x,y} \right)}}.}}}}$

Wherein δ(x, y) is the retardation of said matrix unit (x, y), and θ(x, y) is the fast axis azimuth angle of said matrix unit (x, y). Then the retardation δ(x, y) can be calculated in the range of 0°˜180° as follows:

when √{square root over (V ₁ ²(x, y)+V ₂ ²(x, y))}{square root over (V ₁ ²(x, y)+V ₂ ²(x, y))}≦V ₃(x, y), the δ(x, y)=arcsin(√{square root over (V ₁ ²(x, y)+V ₂ ²(x, y))}{square root over (V ₁ ²(x, y)+V ₂ ²(x, y))}),

when √{square root over (V ₁ ²(x, y)+V ₂ ²(x, y))}{square root over (V ₁ ²(x, y)+V ₂ ²(x, y))}>|V ₃(x, y)|, the δ(x, y)=arccos(V ₃(x, y)),

when √{square root over (V ₁ ²(x, y)+V ₂ ²(x, y))}{square root over (V ₁ ²(x, y)+V ₂ ²(x, y))}<−V ₃(x, y), the δ(x, y)=180°−arcsin(√{square root over (V ₁ ²(x, y)+V ₂ ²(x, y))}{square root over (V ₁ ²(x, y)+V ₂ ²(x, y))}).

The fast axis azimuth angle θ(x, y) can be calculated in the range of −90°˜90° as follows:

when V ₂(x, y)<0&V ₁(x, y)≦0, the

${{\theta \left( {x,y} \right)} = {{\frac{1}{2}{\arctan \left( \frac{V_{1}\left( {x,y} \right)}{V_{2}\left( {x,y} \right)} \right)}} - {90{^\circ}}}},$

when V ₂(x, y)>0, the

${{\theta \left( {x,y} \right)} = {\frac{1}{2}{\arctan \left( \frac{V_{1}\left( {x,y} \right)}{V_{2}\left( {x,y} \right)} \right)}}},$

when V ₂(x, y)<0&V ₁(x, y)>0, the

${\theta \left( {x,y} \right)} = {{\frac{1}{2}{\arctan \left( \frac{V_{1}\left( {x,y} \right)}{V_{2}\left( {x,y} \right)} \right)}} + {90{{^\circ}.}}}$

{circle around (4)} Sequentially change the coordinate values x and y of said matrix unit (x, y) and its corresponding intensity values I₁(x, y), I₂(x, y), I₃(x, y) and I₄(x, y), and then repeat step {circle around (3)}. When the entire matrix units of said measuring sample are calculated, its phase retardation distribution and the fast axis azimuth angle distribution can be obtained.

When the polarization direction of said first analyzer is 135° relative to the fast axis azimuth angle of said quarter wave plate, said computer will perform the steps {circle around (5)}, {circle around (6)}.

{circle around (5)} Said computer processes the intensity values I₁(x, y), I₂(x, y), I₃(x, y) and I₄(x, y) corresponding to said matrix unit (x, y) as follows:

${{V_{1}\left( {x,y} \right)} = {{{\sin \left( {\delta \left( {x,y} \right)} \right)}{\sin \left( {2{\theta \left( {x,y} \right)}} \right)}} = {\frac{2\; {I_{2}\left( {x,y} \right)}}{{I_{2}\left( {x,y} \right)} + {I_{4}\left( {x,y} \right)}} - 1}}},{{V_{2}\left( {x,y} \right)} = {{{\sin \left( {\delta \left( {x,y} \right)} \right)}{\cos \left( {2{\theta \left( {x,y} \right)}} \right)}} = {1 - \frac{2\; {I_{3}\left( {x,y} \right)}}{{I_{2}\left( {x,y} \right)} + {I_{4}\left( {x,y} \right)}}}}},{{V_{3}\left( {x,y} \right)} = {{\cos \left( {\delta \left( {x,y} \right)} \right)} = {\frac{2\; {I_{1}\left( {x,y} \right)}}{{I_{2}\left( {x,y} \right)} + {I_{4}\left( {x,y} \right)}} - 1.}}}$

Then the retardation δ(x, y) can be calculated in the range of 0°˜180° as follows:

when √{square root over (V ₁ ²(x, y)+V ₂ ²(x, y))}{square root over (V ₁ ²(x, y)+V ₂ ²(x, y))}≦V ₃(x, y), the δ(x, y)=arcsin(√{square root over (V ₁ ²(x, y)+V ₂ ²(x, y))}{square root over (V ₁ ²(x, y)+V ₂ ²(x, y))}),

when √{square root over (V ₁ ²(x, y)+V ₂ ²(x, y))}{square root over (V ₁ ²(x, y)+V ₂ ²(x, y))}>|V ₃(x, y)|, the δ(x, y)=arccos(V ₃(x, y)),

when √{square root over (V ₁ ²(x, y)+V ₂ ²(x, y))}{square root over (V ₁ ²(x, y)+V ₂ ²(x, y))}<−V ₃(x, y), the δ(x, y)=180°−arcsin(√{square root over (V ₁ ²(x, y)+V ₂ ²(x, y))}{square root over (V ₁ ²(x, y)+V ₂ ²(x, y))}).

The fast axis azimuth angle θ(x, y) can be calculated in the range of −90°˜90° as follows:

when V ₂(x, y)<0&V ₁(x, y)≦0, the

${{\theta \left( {x,y} \right)} = {{\frac{1}{2}{\arctan \left( \frac{V_{1}\left( {x,y} \right)}{V_{2}\left( {x,y} \right)} \right)}} - {90{^\circ}}}},$

when V ₂(x, y)>0, the

${{\theta \left( {x,y} \right)} = {\frac{1}{2}{\arctan \left( \frac{V_{1}\left( {x,y} \right)}{V_{2}\left( {x,y} \right)} \right)}}},$

when V ₂(x, y)<0&V ₁(x, y)>0, the

${\theta \left( {x,y} \right)} = {{\frac{1}{2}{\arctan \left( \frac{V_{1}\left( {x,y} \right)}{V_{2}\left( {x,y} \right)} \right)}} + {90{{^\circ}.}}}$

{circle around (6)} Sequentially change the coordinate values x and y of said matrix unit (x, y) and its corresponding intensity values I₁(x, y), I₂(x, y), I₃(x, y) and I₄(x, y), and then repeat step {circle around (5)}. When the entire matrix units of said measuring sample are calculated, its phase retardation distribution and the fast axis azimuth angle distribution will be obtained.

Comparing with the previous technology, the technical effect of this invention contains:

1. This invention can measure the phase retardation distribution and fast axis azimuth angle distribution in real time. The light intensity distributions of the four sub-beams are the functions of the phase retardation distribution and fast axis azimuth angle distribution of the measuring sample. The four sub-beams are simultaneously detected by the CCD image sensor and processed at high speed by the computer, thus the phase retardation distribution and fast axis azimuth angle distribution can be obtained in real time.

2. Fluctuation of the initial light intensity will not affect the measured results. The initial light intensity is eliminated during calculation, thus the measured phase retardation distribution and fast axis azimuth angle distribution of the measuring sample is immune to the initial light intensity.

3. The phase retardation and fast axis azimuth angle are of wide measurement range. Using the intensity values of the four sub-beams, the sine and cosine function of the phase retardation and fast axis azimuth angle can be calculated. Utilizing these two functions, the phase retardation can be precisely calculated in the range of 0°˜180° and the fast axis azimuth angle can be precisely calculated in the range of −90°˜90°.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 illustrates the diagram to measure the phase retardation distribution and the fast axis azimuth angle distribution in real time.

FIG. 2 illustrates the diagram of the analyzer array.

DETAILED DESCRIPTION OF THE INVENTION

The invention is further explained in combination with the embodiments and the accompanying diagrams, but the protection scope of the invention should not be limited by it.

EXAMPLE 1

The diagram of measuring the phase retardation distribution and the fast axis azimuth angle distribution in real time is illustrated in FIG. 1. The device of measuring the phase retardation distribution and fast axis azimuth angle distribution in real time is composed of collimating light source 1, circular polarizer 2, diffractive component 4, quarter wave plate 5, analyzer array 6, CCD image sensor 7 and computer 8. Light emitted from the collimating light source 1 successively passes through the circular polarizer 2 and the diffractive component 4 and then is split into four sub-beams. One sub-beam passes through the quarter wave plate 5 and then is analyzed by one analyzer of the analyzer array 6, while the other three sub-beams are directly analyzed by the other three analyzers of the analyzer array 6 without passing through the quarter wave plate 5. The CCD image sensor 7 is connected to the computer 8 through electronics. The faucet for the measuring sample 3 is set between the circular polarizer 2 and the diffractive component 4.

The collimating light source 1 is a He—Ne laser.

The circular polarizer 2 is made of calcite crystal and quartz crystal, whose extinction ratio is better than 10⁻³.

The diffractive component 4 is a Dammann Grating which can split the incident beam into four plus or minus one-order sub-beams, whose light intensity is equal.

The quarter wave plate 5 is a zero-order standard quartz quarter wave plate, which is located in one sub-beam's path generated by diffractive component 4.

The diagram of the analyzer array 6 is illustrated in FIG. 2. It is composed of first analyzer 61, second analyzer 62, third analyzer 63 and fourth analyzer 64. Their polarization direction angles successively increase by 45° and their extinction ratios are all better than 10⁻³. The quarter wave plate 5 is located in the same path with the first analyzer 61. The polarization directions of the first analyzer 61, the second analyzer 62, the third analyzer 63 and the fourth analyzer 64 are respectively 45°, 90°, 135° and 0° relative to the fast axis of the quarter wave plate 5.

The computer 8 is the computer equipped with an image acquisition card.

The method of measuring the phase retardation distribution and fast axis azimuth angle distribution in real time features that it contains the following steps.

{circle around (1)} Insert the measuring sample 3 into the faucet between the circular polarizer 2 and the diffractive component 4, and then adjust the beam to let it perpendicularly pass through the measuring sample 3;

{circle around (2)} Turn on the collimating light source 1, the CCD image sensor 7 and the computer 8.

The CCD image sensor 7 receives an image formed by said four sub-beams and then transmits it to the computer 8. The computer 8 segregates said image into four sub-images. Then the four said sub-images are all pixelated and established a same coordinate system, respectively. The measuring sample 3 is also matrixed and established a coordinate system which is the same with that of the sub-image. The intensity values in said four sub-images corresponding to the measured matrix unit (x, y) of the measuring sample are I₁(x, y), I₂(x, y), I₃(x, y) and I₄(x, y);

{circle around (3)} The computer 8 processes the intensity values I₁(x, y), I₂(x, y), I₃(x, y) and I₄(x, y) corresponding to the matrix unit (x, y) of the measuring sample 3 as follows:

${{V_{1}\left( {x,y} \right)} = {{{\sin \left( {\delta \left( {x,y} \right)} \right)}{\sin \left( {2{\theta \left( {x,y} \right)}} \right)}} = {\frac{2\; {I_{1}\left( {x,y} \right)}}{{I_{1}\left( {x,y} \right)} + {I_{3}\left( {x,y} \right)}} - 1}}},{{V_{2}\left( {x,y} \right)} = {{{\sin \left( {\delta \left( {x,y} \right)} \right)}{\cos \left( {2{\theta \left( {x,y} \right)}} \right)}} = {\frac{2\; {I_{4}\left( {x,y} \right)}}{{I_{1}\left( {x,y} \right)} + {I_{3}\left( {x,y} \right)}} - 1}}},{{V_{3}\left( {x,y} \right)} = {{\cos \left( {\delta \left( {x,y} \right)} \right)} = {1 - {\frac{2\; {I_{2}\left( {x,y} \right)}}{{I_{1}\left( {x,y} \right)} + {I_{3}\left( {x,y} \right)}}.}}}}$

The retardation δ(x, y) of the matrix unit (x, y) can be calculated in the range of 0°˜180°.

When √{square root over (V ₁ ²(x, y)+V ₂ ²(x, y))}{square root over (V ₁ ²(x, y)+V ₂ ²(x, y))}≦V ₃(x, y), the δ(x, y)=arcsin(√{square root over (V ₁ ²(x, y)+V ₂ ²(x, y))}{square root over (V ₁ ²(x, y)+V ₂ ²(x, y))}),

when √{square root over (V ₁ ²(x, y)+V ₂ ²(x, y))}{square root over (V ₁ ²(x, y)+V ₂ ²(x, y))}>|V ₃(x, y), the δ(x, y)=arccos(V ₃(x, y)),

when √{square root over (V ₁ ²(x, y)+V ₂ ²(x, y))}{square root over (V ₁ ²(x, y)+V ₂ ²(x, y))}<−V ₃(x, y), the δ(x, y)=180°−arcsin(√{square root over (V ₁ ²(x, y)+V ₂ ²(x, y))}{square root over (V ₁ ²(x, y)+V ₂ ²(x, y))}).

The fast axis azimuth angle θ(x, y) of the matrix unit (x, y) can be calculated in the range of −90°˜90°.

When V ₂(x, y)<0&V(x, y)≦0, the

${{\theta \left( {x,y} \right)} = {{\frac{1}{2}{\arctan \left( \frac{V_{1}\left( {x,y} \right)}{V_{2}\left( {x,y} \right)} \right)}} - {90{^\circ}}}},$

when V ₂(x, y)>0, the

${{\theta \left( {x,y} \right)} = {\frac{1}{2}{\arctan \left( \frac{V_{1}\left( {x,y} \right)}{V_{2}\left( {x,y} \right)} \right)}}},$

when V ₂(x, y)<0&V ₁(x, y)>0, the

${\theta \left( {x,y} \right)} = {{\frac{1}{2}{\arctan \left( \frac{V_{1}\left( {x,y} \right)}{V_{2}\left( {x,y} \right)} \right)}} + {90{{^\circ}.}}}$

{circle around (4)} Sequentially change the coordinate values x and y of said matrix unit (x, y) and its corresponding intensity values I₁(x, y), I₂(x, y), I₃(x, y) and I₄(x, y), and then repeat step {circle around (3)}. When the entire matrix units of said measuring sample are calculated, its phase retardation distribution and the fast axis azimuth angle distribution will be obtained.

The principle in detail of this invention is explained as follows.

The Jones Vector E of the circularly polarized light emitted from the circular polarizer can be expressed as

$\begin{matrix} {{E = {\frac{E_{0}}{\sqrt{2}}\begin{bmatrix} 1 \\  \end{bmatrix}}},} & (1) \end{matrix}$

wherein E₀ is the amplitude of the circularly polarized light.

The Jones Matrix J_(S) of the matrix unit (x, y) of the measuring sample 3 can be expressed as

$\begin{matrix} {{J_{s} = \begin{bmatrix} {{\cos \frac{\delta \left( {x,y} \right)}{2}} - {\; \sin \frac{\delta \left( {x,y} \right)}{2}\cos \; 2{\theta \left( {x,y} \right)}}} & {{- }\; \sin \frac{\delta \left( {x,y} \right)}{2}\sin \; 2\; {\theta \left( {x,y} \right)}} \\ {{- }\; \sin \frac{\delta \left( {x,y} \right)}{2}\sin \; 2\; {\theta \left( {x,y} \right)}} & {{\cos \frac{\delta \left( {x,y} \right)}{2}} + {\; \sin \frac{\delta \left( {x,y} \right)}{2}\cos \; 2{\theta \left( {x,y} \right)}}} \end{bmatrix}},} & (2) \end{matrix}$

wherein δ(x, y) and θ(x, y) are respectively the phase retardation and the fast axis azimuth angle of the matrix unit (x, y) of the measuring sample 3. The Jones Matrixes J_(P) of the first analyzer 61, the second analyzer 62, the third analyzer 63 and the forth analyzer 64 all can be expressed as

$\begin{matrix} {{J_{P} = \begin{bmatrix} {\cos^{2}\alpha} & {\sin \; {\alpha cos\alpha}} \\ {\sin \; \alpha \; \cos \; \alpha} & {\sin^{2}\alpha} \end{bmatrix}},} & (3) \end{matrix}$

wherein α is the polarization direction angles of the analyzers. The Jones Matrix J_(Q) of the quarter wave plate 5 can be expressed as

$\begin{matrix} {J_{Q} = {\begin{bmatrix} 1 & 0 \\ 0 &  \end{bmatrix}.}} & (4) \end{matrix}$

After directly analyzed by the second analyzer 62, the third analyzer 63 and the fourth analyzer 64 without passing through the quarter wave plate 5, the Jones Vectors E₁(x, y) of three sub-beams all can be expressed as

E ₁(x, y)=J _(P) J _(S) E.   (5)

After passing through the quarter wave plate 5 and then analyzed by the first analyzer 61, the Jones Matrix E₂(x, y) of one sub-beam can be expressed as

E ₂(x, y)=J _(P) J _(Q) J _(S) E.   (6)

When the Jones Matrix E₁(x, y) or E₂(x, y) is multiplied by its conjugate transposed matrix, the intensity values I₁(x, y), I₂(x, y), I₃(x, y), I₄(x, y) of the four sub-beams corresponding to the matrix unit (x, y) of the sample 3 can be expressed as

I ₁(x, y)=I ₀(1−cos(δ(x, y))),   (7)

I ₂(x, y)=I ₀(1−sin(δ(x, y))sin(2θ(x, y))),   (8)

I ₃(x, y)=I ₀(1+sin(δ(x, y))cos(2θ(x, y))),   (9)

I ₄(x, y)=I ₀(1+sin(δ(x, y))sin(2θ(x, y))).   (10)

From equations (7)˜(10), it can be deduced that

$\begin{matrix} {{{V_{1}\left( {x,y} \right)} = {{{\sin \left( {\delta \left( {x,y} \right)} \right)}{\sin \left( {2\; {\theta \left( {x,y} \right)}} \right)}} = {\frac{2{I_{4}\left( {x,y} \right)}}{{I_{2}\left( {x,y} \right)} + {I_{4}\left( {x,y} \right)}} - 1}}},} & (11) \\ {{{V_{2}\left( {x,y} \right)} = {{{\sin \left( {\delta \left( {x,y} \right)} \right)}{\cos \left( {2\; {\theta \left( {x,y} \right)}} \right)}} = {\frac{2{I_{3}\left( {x,y} \right)}}{{I_{2}\left( {x,y} \right)} + {I_{4}\left( {x,y} \right)}} - 1}}},} & (12) \\ {{V_{3}\left( {x,y} \right)} = {{\cos \left( {\delta \left( {x,y} \right)} \right)} = {1 - {\frac{2{I_{1}\left( {x,y} \right)}}{{I_{2}\left( {x,y} \right)} + {I_{4}\left( {x,y} \right)}}.}}}} & (13) \end{matrix}$ Then when √{square root over (V ₁ ²(x, y)+V ₂ ²(x, y))}{square root over (V ₁ ²(x, y)+V ₂ ²(x, y))}≦V ₃(x, y), the δ(x, y)=arcsin(√{square root over (V ₁ ²(x, y)+V ₂ ²(x, y))}{square root over (V ₁ ²(x, y)+V ₂ ²(x, y))}),   (14)

when √{square root over (V ₁ ²(x, y)+V ₂ ²(x, y))}{square root over (V ₁ ²(x, y)+V ₂ ²(x, y))}>|V ₃(x, y)|, the δ(x, y)=arccos(V ₃(x, y)),   (15)

when √{square root over (V ₁ ²(x, y)+V ₂ ²(x, y))}{square root over (V ₁ ²(x, y)+V ₂ ²(x, y))}<−V ₃(x, y), the δ(x, y)=180°−arcsin(√{square root over (V ₁ ²(x, y)+V ₂ ²(x, y))}{square root over (V ₁ ²(x, y)+V ₂ ²(x, y))}).   (16)

When V ₂(x, y)<0&V ₁(x, y)≦0, the

$\begin{matrix} {{{\theta \left( {x,y} \right)} = {{\frac{1}{2}{\arctan \left( \frac{V_{1}\left( {x,y} \right)}{V_{2}\left( {x,y} \right)} \right)}} - {90{^\circ}}}},} & (17) \end{matrix}$

when V ₂(x, y)>0, the

$\begin{matrix} {{{\theta \left( {x,y} \right)} = {\frac{1}{2}{\arctan \left( \frac{V_{1}\left( {x,y} \right)}{V_{2}\left( {x,y} \right)} \right)}}},} & (18) \end{matrix}$

when V ₂(x, y)<0&V ₁(x, y)>0, the

$\begin{matrix} {{\theta \left( {x,y} \right)} = {{\frac{1}{2}{\arctan \left( \frac{V_{1}\left( {x,y} \right)}{V_{2}\left( {x,y} \right)} \right)}} + {90{{^\circ}.}}}} & (19) \end{matrix}$

y using equations (14)˜(19), the δ(x, y) distribution can be measured in the range of 0°˜180° and θ(x, y) distribution can be measured in the range of −90°˜90°.

EXAMPLE 2

The difference between embodiment 2 and embodiment 1 is that the polarization directions of the first analyzer 61, the second analyzer 62, the third analyzer 63 and the fourth analyzer 64 are 135°, 0°, 45° and 90° relative to the fast axis of the quarter wave plate 5. And the corresponding data processing steps are also different.

{circle around (3)} The computer 8 processes the intensity values I₁(x, y), I₂(x, y), I₃(x, y) and I₄(x, y) corresponding to the matrix unit (x, y) of the measuring sample 3 as follows:

$\begin{matrix} {{{V_{1}\left( {x,y} \right)} = {{{\sin \left( {\delta \left( {x,y} \right)} \right)}{\sin \left( {2\; {\theta \left( {x,y} \right)}} \right)}} = {\frac{2{I_{2}\left( {x,y} \right)}}{{I_{2}\left( {x,y} \right)} + {I_{4}\left( {x,y} \right)}} - 1}}},} \\ {{{V_{2}\left( {x,y} \right)} = {{{\sin \left( {\delta \left( {x,y} \right)} \right)}{\cos \left( {2\; {\theta \left( {x,y} \right)}} \right)}} = {1 - \frac{2{I_{3}\left( {x,y} \right)}}{{I_{2}\left( {x,y} \right)} + {I_{4}\left( {x,y} \right)}}}}},} \\ {{V_{3}\left( {x,y} \right)} = {{\cos \left( {\delta \left( {x,y} \right)} \right)} = {\frac{2{I_{1}\left( {x,y} \right)}}{{I_{2}\left( {x,y} \right)} + {I_{4}\left( {x,y} \right)}} - 1.}}} \end{matrix}$

The retardation δ(x, y) of the matrix unit (x, y) can be calculated in the range of 0˜180°.

When √{square root over (V ₁ ²(x, y)+V ₂ ²(x, y))}{square root over (V ₁ ²(x, y)+V ₂ ²(x, y))}≦V ₃(x, y), the δ(x, y)=arcsin(√{square root over (V ₁ ²(x, y)+V ₂ ²(x, y))}{square root over (V ₁ ²(x, y)+V ₂ ²(x, y))}),

when √{square root over (V ₁ ²(x, y)+V ₂ ²(x, y))}{square root over (V ₁ ²(x, y)+V ₂ ²(x, y))}>|V ₃(x, y)|, the δ(x, y)=arccos(V ₃(x, y)),

when √{square root over (V ₁ ²(x, y)+V ₂ ²(x, y))}{square root over (V ₁ ²(x, y)+V ₂ ²(x, y))}<−V ₃(x, y), the δ(x, y)=180°−arcsin(√{square root over (V ₁ ²(x, y)+V ₂ ²(x, y))}{square root over (V ₁ ²(x, y)+V ₂ ²(x, y))}).

The fast axis azimuth angle θ(x, y) of the matrix unit (x, y) can be calculated in the range of −90°˜90°.

When V ₂(x, y)<0&V ₁(x, y)≦0, the

${{\theta \left( {x,y} \right)} = {{\frac{1}{2}{\arctan \left( \frac{V_{1}\left( {x,y} \right)}{V_{2}\left( {x,y} \right)} \right)}} - {90{^\circ}}}},$

when V ₂(x, y)>0, the

${{\theta \left( {x,y} \right)} = {\frac{1}{2}{\arctan \left( \frac{V_{1}\left( {x,y} \right)}{V_{2}\left( {x,y} \right)} \right)}}},$

when V ₂(x, y)<0&V ₁(x, y)>0, the

${\theta \left( {x,y} \right)} = {{\frac{1}{2}{\arctan \left( \frac{V_{1}\left( {x,y} \right)}{V_{2}\left( {x,y} \right)} \right)}} + {90{{^\circ}.}}}$

{circle around (4)} Sequentially change the coordinate values x and y of said matrix unit (x, y) and its corresponding intensity values I₁(x, y), I₂(x, y), I₃(x, y) and I₄(x, y), and then repeat step {circle around (3)}. When the entire matrix units of said measuring sample are calculated, its phase retardation distribution and the fast axis azimuth angle distribution will be obtained.

The principle of this invention is explained as follows.

The Jones Vector E of the circularly polarized light emitted from the circular polarizer can be expressed as

$\begin{matrix} {{E = {\frac{E_{0}}{\sqrt{2}}\begin{bmatrix} 1 \\  \end{bmatrix}}},} & (20) \end{matrix}$

wherein E₀ is the amplitude of the circularly polarized light. The Jones Matrix J_(S) of the matrix unit (x, y) of the measuring sample 3 can be expressed as

$\begin{matrix} {{J_{s} = \begin{bmatrix} {{\cos \frac{\delta \left( {x,y} \right)}{2}} - {i\; \sin \frac{\delta \left( {x,y} \right)}{2}\cos \; 2{\theta \left( {x,y} \right)}}} & {{- i}\; \sin \frac{\delta \left( {x,y} \right)}{2}\sin \; 2{\theta \left( {x,y} \right)}} \\ {{- i}\; \sin \frac{\delta \left( {x,y} \right)}{2}\sin \; 2{\theta \left( {x,y} \right)}} & {{\cos \frac{\delta \left( {x,y} \right)}{2}} + {i\; \sin \frac{\delta \left( {x,y} \right)}{2}\cos \; 2{\theta \left( {x,y} \right)}}} \end{bmatrix}},} & (21) \end{matrix}$

wherein δ(x, y) and θ(x, y) are respectively the phase retardation and the fast axis azimuth angle of the matrix unit (x, y) of the sample 3. The Jones Matrixes J_(P) of the first analyzer 61, the second analyzer 62, the third analyzer 63 and the forth analyzer 64 all can be expressed as

$\begin{matrix} {{J_{P} = \begin{bmatrix} {\cos^{2}\alpha} & {\sin \; {\alpha cos\alpha}} \\ {\sin \; {\alpha cos\alpha}} & {\sin^{2}\alpha} \end{bmatrix}},} & (22) \end{matrix}$

wherein α is the polarization direction angles of the analyzers. The Jones Matrix J_(Q) of the quarter wave plate 5 can be expressed as

$J_{Q} = {\begin{bmatrix} 1 & 0 \\ 0 & i \end{bmatrix}.}$

After directly analyzed by the second analyzer 62, the third analyzer 63 and the fourth analyzer 64 without passing through the quarter wave plate 5, the Jones Vectors E₁(x, y) of three sub-beams all can be expressed as

E ₁(x, y)=J _(P) J _(S) E.   (24)

After passing through the quarter wave plate 5 and then analyzed by the first analyzer 61, the Jones Matrix E₂(x, y) of one sub-beam can be expressed as

E ₂(x, y)=J _(P) J _(Q) J _(S) E.   (25)

When the Jones Matrix E₁(x, y) or E₂(x, y) is multiplied by its conjugate transposed matrix, the intensity values I₁(x, y), I₂(x, y), I₃(x, y), I₄(x, y) of the four sub-beams corresponding to the matrix unit (x, y) of the measuring sample 3 can be expressed as

I ₁(x, y)=I ₀(1+cos(δ(x, y))),   (26)

I ₂(x, y)=I ₀(1+sin(δ(x, y))sin(2θ(x, y))),   (27)

I ₃(x, y)=I ₀(1−sin(δ(x, y))cos(2θ(x, y))),   (28)

I ₄(x, y)=I ₀(1−sin(δ(x, y))sin(2θ(x, y))).   (29)

From equations (26)˜(29), it can be deduced that

$\begin{matrix} \begin{matrix} {{V_{1}\left( {x,y} \right)} = {{\sin \left( {\delta \left( {x,y} \right)} \right)}{\sin \left( {2{\theta \left( {x,y} \right)}} \right)}}} \\ {{= {\frac{2{I_{2}\left( {x,y} \right)}}{{I_{2}\left( {x,y} \right)} + {I_{4}\left( {x,y} \right)}} - 1}},} \end{matrix} & (30) \\ \begin{matrix} {{V_{2}\left( {x,y} \right)} = {{\sin \left( {\delta \left( {x,y} \right)} \right)}{\cos \left( {2{\theta \left( {x,y} \right)}} \right)}}} \\ {{= {1 - \frac{2{I_{3}\left( {x,y} \right)}}{{I_{2}\left( {x,y} \right)} + {I_{4}\left( {x,y} \right)}}}},} \end{matrix} & (31) \\ \begin{matrix} {{V_{3}\left( {x,y} \right)} = {\cos \left( {\delta \left( {x,y} \right)} \right)}} \\ {= {\frac{2{I_{1}\left( {x,y} \right)}}{{I_{2}\left( {x,y} \right)} + {I_{4}\left( {x,y} \right)}} - 1.}} \end{matrix} & (32) \end{matrix}$ Then when √{square root over (V ₁ ²(x, y)+V ₂ ²(x, y))}{square root over (V ₁ ²(x, y)+V ₂ ²(x, y))}≦V ₃(x, y), the δ(x, y)=arcsin(√{square root over (V ₁ ²(x, y)+V ₂ ²(x, y))}{square root over (V ₁ ²(x, y)+V ₂ ²(x, y))}),   (33)

when √{square root over (V ₁ ²(x, y)+V ₂ ²(x, y))}{square root over (V ₁ ²(x, y)+V ₂ ²(x, y))}>|V ₃(x, y)|, the δ(x, y)=arccos(V ₃(x, y)),   (34)

when √{square root over (V ₁ ²(x, y)+V ₂ ²(x, y))}{square root over (V ₁ ²(x, y)+V ₂ ²(x, y))}<−V ₃(x, y), the δ(x, y)=180°−arcsin(√{square root over (V ₁ ²(x, y)+V ₂ ²(x, y))}{square root over (V ₁ ²(x, y)+V ₂ ²(x, y))}).   (35)

When V ₂(x, y)<0&V ₁(x, y)≦0, the

$\begin{matrix} {{{\theta \left( {x,y} \right)} = {{\frac{1}{2}{\arctan \left( \frac{V_{1}\left( {x,y} \right)}{V_{2}\left( {x,y} \right)} \right)}} - {90{^\circ}}}},} & (36) \end{matrix}$

when V ₂(x, y)>0, the

$\begin{matrix} {{{\theta \left( {x,y} \right)} = {\frac{1}{2}{\arctan \left( \frac{V_{1}\left( {x,y} \right)}{V_{2}\left( {x,y} \right)} \right)}}},} & (37) \end{matrix}$

when V ₂(x, y)<0&V ₁(x, y)>0, the

$\begin{matrix} {{\theta \left( {x,y} \right)} = {{\frac{1}{2}{\arctan \left( \frac{V_{1}\left( {x,y} \right)}{V_{2}\left( {x,y} \right)} \right)}} + {90{{^\circ}.}}}} & (38) \end{matrix}$

By using equations (33)˜(38), the δ(x, y) distribution can be measured in the range of 0°˜180° and θ(x, y) distribution can be measured in the range of −90°˜90°.

Experiment results show that this invention can measure the phase retardation distribution and fast axis azimuth angle distribution in real time. The measured results are immune to the fluctuation of the initial intensity. Meanwhile this invention is of wide measurement range. 

We claim:
 1. A device for measuring the phase retardation distribution and fast axis azimuth angle distribution in real time, comprising a collimating light source, a circular polarizer, a diffractive beam-splitting component, a quarter wave plate, an analyzer array, a CCD image sensor, and a computer equipped with an image acquisition card, wherein the analyzer array is composed of four analyzers whose polarization direction angles successively increase by 45°, they are respectively named as a first analyzer, a second analyzer, a third analyzer, and a fourth analyzer; the quarter wave plate is located in the same path with the first analyzer, and the angle between a fast axis of the quarter wave plate and transmission direction of the first analyzer is 45° or 135°; light emitted from the collimating light source passes through the circular polarizer and the diffractive component, and is then split into four sub-beams; one sub-beam passes through the quarter wave plate and then is analyzed by the first analyzer, the other three sub-beams are directly analyzed by the second analyzer, the third analyzer, and the fourth analyzer, respectively; an output port of the image sensor is connected to an input port of the computer; and a faucet for the measuring sample is set between the circular polarizer and the diffractive component.
 2. The device for measuring the phase retardation distribution and fast axis azimuth angle distribution in real time of claim 1, wherein the collimating light source is a He—Ne laser, solid state laser, or semiconductor laser.
 3. The device for measuring the phase retardation distribution and fast axis azimuth angle distribution in real time of claim 1, wherein the circular polarizer is made of birefringent crystal, birefringent film, or micro optical element.
 4. The device for measuring the phase retardation distribution and fast axis azimuth angle distribution in real time of claim 1, wherein the diffractive component is Quadrature Amplitude grating, Quadrature Phase grating or Dammann grating, which splits the incident beam into four sub-beams with the same light intensity.
 5. The device for measuring the phase retardation distribution and fast axis azimuth angle distribution in real time of claim 1, wherein the quarter wave plate is a crystal wave plate, prismatic wave plate, film wave plate, or composite wave plate.
 6. The device for measuring the phase retardation distribution and fast axis azimuth angle distribution in real time of claim 1, wherein the first analyzer, second analyzer, third analyzer, and fourth analyzer are all the polarizer whose extinction ratio is better than 10⁻³.
 7. The method for measuring the phase retardation distribution and fast axis azimuth angle distribution in real time of claim 1, comprising inserting a measuring sample into the faucet between said circular polarizer and the diffractive component, adjusting the light beam to let it perpendicularly pass through the measuring sample, forming an image by the four sub-beams, receiving the image by the CCD image sensor and transmitting the image to the computer, segregating the image into four sub-images in the computer, and processing the four sub-images in the computer to obtain a phase retardation distribution and fast axis azimuth angle distribution of the measuring sample, wherein the four sub-images are all pixelated and established a same coordinate system, respectively; the measuring sample is also matrixed and established a coordinate system which is the same with the coordinate system of the sub-image; intensity values in said four sub-images corresponding to the measured matrix unit (x, y) of the measuring sample are I₁(x, y), I₂(x, y), I₃(x, y) and I₄(x, y); calculating said intensity values to obtain the phase retardation and fast axis azimuth angle of said matrix unit (x, y) of the measuring sample.
 8. The measurement method of claim 7, wherein the computer processes said sub-images by when the polarization direction of the first analyzer is 45° relative to the fast axis of the quarter wave plate, performing steps of (3) and (4) by the computer: (3) processing the intensity values I₁(x, y), I₂(x, y), I₃(x, y) and I₄(x, y) corresponding to said matrix unit (x, y) as follows: $\begin{matrix} {{V_{1}\left( {x,y} \right)} = {{\sin \left( {\delta \left( {x,y} \right)} \right)}{\sin \left( {2{\theta \left( {x,y} \right)}} \right)}}} \\ {{= {\frac{2{I_{1}\left( {x,y} \right)}}{{I_{1}\left( {x,y} \right)} + {I_{3}\left( {x,y} \right)}} - 1}},} \end{matrix}$ $\begin{matrix} {{V_{2}\left( {x,y} \right)} = {{\sin \left( {\delta \left( {x,y} \right)} \right)}{\cos \left( {2{\theta \left( {x,y} \right)}} \right)}}} \\ {{= {\frac{2{I_{4}\left( {x,y} \right)}}{{I_{1}\left( {x,y} \right)} + {I_{3}\left( {x,y} \right)}} - 1}},} \end{matrix}$ $\begin{matrix} {{V_{3}\left( {x,y} \right)} = {\cos \left( {\delta \left( {x,y} \right)} \right)}} \\ {= {1 - {\frac{2{I_{2}\left( {x,y} \right)}}{{I_{1}\left( {x,y} \right)} + {I_{3}\left( {x,y} \right)}}.}}} \end{matrix}$ wherein δ(x, y) is a retardation of the matrix unit (x, y), and θ(x, y) is a fast axis azimuth angle of the matrix unit (x, y); calculating the retardation δ(x, y) in a range of 0°˜180° as follows: when √{square root over (V ₁ ²(x, y)+V ₂ ²(x, y))}{square root over (V ₁ ²(x, y)+V ₂ ²(x, y))}≦V ₃(x, y), the δ(x, y)=arcsin(√{square root over (V ₁ ²(x, y)+V ₂ ²(x, y))}{square root over (V ₁ ²(x, y)+V ₂ ²(x, y))}), when √{square root over (V ₁ ²(x, y)+V ₂ ²(x, y))}{square root over (V ₁ ²(x, y)+V ₂ ²(x, y))}>|V ₃(x, y)|, the δ(x, y)=arccos(V ₃(x, y)), when √{square root over (V ₁ ²(x, y)+V ₂ ²(x, y))}{square root over (V ₁ ²(x, y)+V ₂ ²(x, y))}<−V ₃(x, y), the δ(x, y)=180°−arcsin(√{square root over (V ₁ ²(x, y)+V ₂ ²(x, y))}{square root over (V ₁ ²(x, y)+V ₂ ²(x, y))}). calculating the fast axis azimuth angle θ(x, y) in a range of −90°˜90° as follows: when V ₂(x, y)<0&V ₁(x, y)≦0, the ${{\theta \left( {x,y} \right)} = {{\frac{1}{2}{\arctan \left( \frac{V_{1}\left( {x,y} \right)}{V_{2}\left( {x,y} \right)} \right)}} - {90{^\circ}}}},$ when V ₂(x, y)>0, the ${{\theta \left( {x,y} \right)} = {\frac{1}{2}{\arctan \left( \frac{V_{1}\left( {x,y} \right)}{V_{2}\left( {x,y} \right)} \right)}}},$ when V ₂(x, y)<0&V ₁(x, y)>0, the ${{\theta \left( {x,y} \right)} = {{\frac{1}{2}{\arctan \left( \frac{V_{1}\left( {x,y} \right)}{V_{2}\left( {x,y} \right)} \right)}} + {90{^\circ}}}},$ and (4) sequentially changing the coordinate values x and y of the matrix unit (x, y) and its corresponding intensities I₁(x, y), I₂(x, y), I₃(x, y) and I₄(x, y), and repeating step (3); and when the entire matrix units of the measuring sample are calculated, its phase retardation distribution and the fast axis azimuth angle distribution are obtained; when the polarization direction of the first analyzer is 135° relative to the fast axis azimuth angle of the quarter wave plate, performing steps (5) and (6) by the computer: (5) processing the intensity values I₁(x, y), I₂(x, y), I₃(x, y) and I₄(x, y) corresponding to said matrix unit (x, y) as follows: $\begin{matrix} {{V_{1}\left( {x,y} \right)} = {{\sin \left( {\delta \left( {x,y} \right)} \right)}{\sin \left( {2{\theta \left( {x,y} \right)}} \right)}}} \\ {{= {\frac{2{I_{2}\left( {x,y} \right)}}{{I_{2}\left( {x,y} \right)} + {I_{4}\left( {x,y} \right)}} - 1}},} \end{matrix}$ $\begin{matrix} {{V_{2}\left( {x,y} \right)} = {{\sin \left( {\delta \left( {x,y} \right)} \right)}{\cos \left( {2{\theta \left( {x,y} \right)}} \right)}}} \\ {{= {1 - \frac{2{I_{3}\left( {x,y} \right)}}{{I_{2}\left( {x,y} \right)} + {I_{4}\left( {x,y} \right)}}}},} \end{matrix}$ $\begin{matrix} {{V_{3}\left( {x,y} \right)} = {\cos \left( {\delta \left( {x,y} \right)} \right)}} \\ {{= {\frac{2{I_{1}\left( {x,y} \right)}}{{I_{2}\left( {x,y} \right)} + {I_{4}\left( {x,y} \right)}} - 1}},} \end{matrix}$ and calculating the retardation δ(x, y) in a range of 0°˜180° as follows: when √{square root over (V ₁ ²(x, y)+V ₂ ²(x, y))}{square root over (V ₁ ²(x, y)+V ₂ ²(x, y))}≧V ₃(x, y), the δ(x, y)=arcsin(√{square root over (V ₁ ²(x, y)+V ₂ ²(x, y))}{square root over (V ₁ ²(x, y)+V ₂ ²(x, y))}), when √{square root over (V ₁ ²(x, y)+V ₂ ²(x, y))}{square root over (V ₁ ²(x, y)+V ₂ ²(x, y))}>|V ₃(x, y)|, the δ(x, y)=arccos(V ₃(x, y)), when √{square root over (V ₁ ²(x, y)+V ₂ ²(x, y))}{square root over (V ₁ ²(x, y)+V ₂ ²(x, y))}<−V ₃(x, y), the δ(x, y)=180°−arcsin(√{square root over (V ₁ ²(x, y)+V ₂ ²(x, y))}{square root over (V ₁ ²(x, y)+V ₂ ²(x, y))}); and calculating the fast axis azimuth angle θ(x, y) in a range of −90°˜90° as follows: when V ₂(x, y)<0&V ₁(x, y)≦0, the ${{\theta \left( {x,y} \right)} = {{\frac{1}{2}{\arctan \left( \frac{V_{1}\left( {x,y} \right)}{V_{2}\left( {x,y} \right)} \right)}} - {90{^\circ}}}},$ when V ₂(x, y)>0, the ${{\theta \left( {x,y} \right)} = {\frac{1}{2}{\arctan \left( \frac{V_{1}\left( {x,y} \right)}{V_{2}\left( {x,y} \right)} \right)}}},$ when V ₂(x, y)<0&V ₁(x, y)>0, the ${{\theta \left( {x,y} \right)} = {{\frac{1}{2}{\arctan \left( \frac{V_{1}\left( {x,y} \right)}{V_{2}\left( {x,y} \right)} \right)}} + {90{^\circ}}}};$ and (6) sequentially changing the coordinate values x and y of the matrix unit (x, y) and its corresponding intensities I₁(x, y), I₂(x, y), I₃(x, y) and I₄(x, y), and repeating step (5), and when the entire matrix units of the measuring sample are calculated, its phase retardation distribution and the fast axis azimuth angle distribution are obtained. 